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Abstract 

The subfilter-scale (SFS) physics of regularization models are in- 
vestigated to understand the regularizations' performance as SFS 
models. Suppression of spectrally local SFS interactions and conser- 
vation of small-scale circulation in the Lagrangian-averaged Navier- 
Stokes model (LANS— a) is found to lead to the formation of rigid 
bodies. These contaminate the superfilter-scale energy spectrum 
with a scaling that approaches k^^ as the SFS spectra is resolved. 
The Clark— a and Leray— a models, truncations of LANS— a, do 
not conserve small-scale circulation and do not develop rigid bodies. 
LANS— a, however, is closest to Navier-Stokes in intermittency prop- 
erties. AH three models are found to be stable at high Reynolds num- 
ber. Differences between and norm models are clarified. For 
magnetohydrodynamics (MHD), the presence of the Lorentz force as 
a source (or sink) for circulation and as a facilitator of both spectrally 
nonlocal large to small scale interactions as well as local SFS interac- 
tions prevents the formation of rigid bodies in Lagrangian-averaged 
MHD (LAMHD-a). LAMHD-q performs weU as a predictor of 
superfilter-scale energy spectra and of intermittent current sheets 
at high Reynolds numbers. It may prove generally applicable as a 
MHD-LES. 
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1 Introduction 



Computing solutions to the Navier-Stokes equations at realistic values of 
the Reynolds number (Re = UL^/v, with characteristic velocity, C/, and 
length-scale, Lq and viscosity, v) for most geophysical, astrophysical, and 
many engineering applications is technologically infeasible. This is because 
the range of dynamically important length (and time) scales is quite large: 
from the large scale, Lq, defined by the problem down to the scales of vis- 
cous dissipation, li, LqRc^^^^ . One approach is to simply cut off the 
smallest scales to arrive at a problem small enough for modern computa- 
tional limits. A low-bandpass filter, L : / ^ /, replaces the velocity, v, and 
pressure, P, with smoother, resolvable fields, v , P. Their time evolution is 
governed by the filtered Navier-Stokes equations, 

9tv -I- a; X V = - VP + j/V^v - V • f 

V-v = 0, (1) 

where a; = V x v is the vorticity and f is the Reynolds subfilter-scale (SFS) 
stress tensor, 

fij = Wvj - ViVj . (2) 

The modeling of the unresolved stress, Eq. ([2]), is the main challenge of 
conducting such large eddy simulations (LES, see [12] for a recent review). 

Regularization modeling (of the SFS stress tensor) for Navier-Stokes 
[Tni[2i[Tl[nilHl[Il|31lini[Sn], magnetohydrodynamlcs (MHD) [33], 
Boussinesq convection [53] , and inviscid cases [37] promises several advan- 
tages. For Navier-Stokes, only weak, possibly non-unique solutions have 
been rigorously proven to exist, and this can impact the possibility of 
achieving a direct numerical solution (DNS), e.g., with Fourier methods 
|26j . This is because to prove, even in the linear case, the convergence of a 
numerical solver for a partial differential equation (PDF), one assumes the 
PDF is well-posed (i.e., the existence of unique solutions that depend con- 
tinuously on the data) [55] • A regularization has strong, unique solutions 
(it is well-posed) and the concerns of |26| do not prevent us from achieving a 
DNS of the model equations, even with Fourier methods. It is worth noting 
that numerical convergence of the DNS of the regularization model equa- 
tions implies a grid-independent model of the Navier-Stokes equations. Ad- 
ditionally, the dissipative term will be unmodified (unlike many LES); the 
Reynolds number remains well defined (e.g., a; x v/i^V^v ~ ULq/v = Re). 
This is opposed to the usual approach of modeling the behavior of the flow 
in the limit of very large Re. Thus, the models may be more applicable 
to intermittent phenomena where the length of the inertial range can be 
important [23] . Since the models do not introduce the effect of the small 
scales in an ad hoc fashion but rather preserve the mathematical properties 
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of the underlying equations, their apphcation can further our understand- 
ing of turbulence and turbulence modeling. The methods are also more 
easily generalized to other problems (e.g., coupling to a magnetic field). 

In this paper we address two separate questions. One is the question of 
the practical applicability of regularization models as SFS models. When 
addressing this question, the filter width, a, will be placed in the inertial 
range and the grid spacing, A, will be just small enough to achieve a DNS 
of the regularization. Our aim is then to determine how well the model's 
DNS reproduces a "DNS" of Navier-Stokes compared at scales larger than 
a. Our second question is "How do the models work?" To answer this 
question, we take an approach that is the antithesis of a practical LES: 
we choose the filter width, a, to be a large fraction of the computation 
domain (and, thus, a large multiple of A). Such a calculation is not a 
LES-there is no superfilter-scale inertial range to compare with Navier- 
Stokes. Instead, we make a new type of study to understand the new 
SFS physics introduced by the regularization. The differences in physics 
between the SFS model and Navier-Stokes (or MHD) is what allows the 
model to reduce computational cost when it is employed as a SFS model. 
Understanding of how the models work (or fail) can guide the development 
of new models as we will show. 

2 Navier-Stokes 

2.1 LANS— a and rigid body formation 

The first model we consider is the Lagrangian-averaged Navier-Stokes (LANS) 
a— model [IUl[2n]. It is derived by Lagrangian averaging fluid motions fol- 
lowed by application of Taylor's frozen-in turbulence approximation as the 
model's one and only closure: fluctuations about the Lagrangian mean 
smaller than a are swept along by the large-scale flow and are not allowed 
to interact with one another [28]. The model is attractive as it retains the 
Hamiltonian structure of Navier-Stokes, preserves Kelvin's theorem (con- 
serves small-scale circulation in the absence of dissipation), and conserves 
both total energy and helicity (the correlation between the velocity, v, and 
its curl, the vorticity a; = V x v) [29]. These properties are conserved 
in the norm instead of the usual norm. This is essential when 
interpreting results of the model as, for example, quantities involving the 
square velocity, |vp, must now be replaced with the dot product v-v where 
V is the filtered velocity. Physically, due to the frozen-in approximation, 
the model retains spectrally nonlocal interactions (important at finite Re 
[43l [2] ) between the superfilter and subfilter scales while the fiux of energy 
in subfilter scales is reduced by the limit on local small-scale to small-scale 
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interactions [Till^ . 

The LANS— a model is given by, 

9tv + u; X V = — Vtt + vV^v 

V • V = V • V = . (3) 

From the identity, a; x v = v • Vv + (Vv)"^ • v — V(v • v), we see that it 
differs from Navier-Stokes both in advection by the smoothed velocity and 
the addition of a second nonlinear term associated with the conservation 
of circulation. Traditionally, LANS— a is used with an inverse Helmholtz 
operator as the filter: Vi = {\ — a'^djj)~^Vi. In this case, LANS— a can be 
written as a LES, Eq. ([T]), with 

= (1 - a^djjy^a^idmVidrnVj + dynVtdjVm - diVmdjVm) ■ (4) 

The model allows for a reduction in resolution without changing (or sup- 
plementing) the dissipative term by instead altering the SFS scaling prop- 
erties. Near wavenumber, k = 27r/a, the energy spectrum is predicted 
to transition from with /3 = —5/3 at larger scales to /? = — 1 at smaller 
scales [2]. Consequently, dissipation (f2(Q)(fc) — k'^E(^a){k)) goes as k^ 
instead of k^f^ and the same amount of power is dissipated in fewer scales. 
The change in spectral scaling also predicts a resolution requirement in 
degrees of freedom, do/, for LANS [M] . 

do/a - a-'^Re''^ , (5) 

which has been confirmed in numerical experiments [22]. Once do fa has 
been resolved, further resolution yields no change in the numerical solution: 
LANS— a is a grid-independent SFS model. When compared with the dof 
for Navier-Stokes, 

dofNs - Re^'^ , (6) 

we see that LANS— a should improve as a SFS model for larger Re. This 
was an encouraging prediction as LANS— a compared well with dynamic 
eddy viscosity |46| and dynamic mixed (similarity) eddy viscosity |18| at 
moderate Re. 

We have, however, found that LANS— a develops a problem at large 
Re: it develops a positive-exponent power- law bump in its small-scale en- 
ergy spectrum and a contamination of superfilter-scale spectral properties 
P^ . To investigate the SFS physics responsible for this, we employ (as the 
antithesis of a LES) a filter 1/3 the size of our 256^ computational cube 
in a pseudo-spectral calculation [3D1 [H] with a Taylor-Green (TG) forcing 
[5T| and Re ~ 8000. As shown in Fig. [1] the observed scaling law is A:"*"^. 
This was shown to be associated with the formation in the flow of pas- 
sively swept regions, called rigid bodies [22 . These form as a consequence 
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Figure 1: Energy spectrum for LANS— q (solid line) with filter scale a = 
27r/3 (vertical dashed line). The SFS power law is well approximated by 
k^^ . Insert shows thresholded cubed velocity increment |(5i;|| (27r/10)p < 
10~^ in black. These regions do not contribute to the turbulent cascade of 
energy to smaller scales and are identified with rigid bodies. A spectrum 
of only the white regions (dashed line) is consistent with the predicted 
scaling outside rigid bodies. 
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of disallowing sub-a-scale fluctuations to interact with each other in the 
closure approximation. A rigid body cannot support longitudinal velocity 
increments: Sv^^{l) = [v(x) — v(x + l)] - l/l = 0. This dimensionally predicts 
a scaling relation, Sv ~ l'^, and, with v ^ {1 + a'^k'^)v ~ a^k'^v for I -C a, 
an energy spectrum of 

Eo,{k) - vvk-^ - k^ (7) 

which is compatible with the observed SFS energy spectrum. Inside rigid 
bodies there can be no turbulent cascade of energy to smaller scales (since 
there are no internal degrees of freedom) . From the Karman-Howarth the- 
orem, we should then expect to be able to detect rigid bodies by visualizing 
the cubed velocity increments (which are proportional to the energy flux). 
The regions which correspond to negligible flux are shown as black in the 
inset of Fig. [TJ Filtering these regions out, allows us to obtain a (convolved) 
energy spectra for the remaining white portion of the flow. This spectrum 
is shown as a dashed line in Fig. [T]and has a negative spectral slope close to 
the predicted k~^ spectrum. The resulting picture of the model's behavior 
is to produce two spatially separate scalings. The white portions of the flow 
possess the predicted LANS— a scaling and are responsible for the obser- 
vation of the predicted do fa- The black portions are rigid bodies whose k^ 
energy spectrum dominates over k~^ for large k and are responsible for the 
observed spectral contamination. Because of this spectral contamination, 
there is no need for further tests of LANS— a against newer models such 
as the Variational Multiscale Method [3^. Note, however, that suitable 
spectra can be obtained with a modified viscous length scale (LANS— a/? 
[33]). 

2.2 The influence of circulation on rigid bodies 

The formation of rigid bodies in LANS— a limits the reduction of numerical 
dof saved compared to Navier-Stokes to a factor of 1/12 regardless of Re 
[22] . It is desirable, then, to alter the model in such a way to prevent the 
formation of rigid bodies. Truncation of the SFS stress tensor, Eq. (|4]), to 
the first term results in the Clark— a model [S], 

Stv + (1 - ]^a^y^){u3 X v) - [(y^u;) x v + (X) x (V^v)] = 

V-v = V-v = 0, (8) 

and to the first two terms results in the Leray— a model [TT l 118 ) . 

9tv + V • Vv = -Vp + i^V^v 

V • V = V ■ V = . (9) 
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Both these models are regularizations and conserve the total energy of the 
flow [8l[T7l[T8]. They do not, however, conserve the helicity nor the small- 
scale circulation. Considering the rotational properties of a rigid body (in 
the absence of viscous friction), these models' circulation properties may be 
incompatible with rigid body formation. Indeed, while LANS— a exhibits 
a positive-exponent power law in this case (a = 27r/13, Re « 3300, TG 
forcing), both Clark— a and Leray— a are free from this signature of rigid 
body formation (Fig. However, LANS— a's intermittency properties 

are more similar to Navier-Stokes than the other two models (Fig. [2}3). 

2.3 Computability and interpretation of norm reg- 
ularizations 

It has been reported recently that Clark-a and LANS-a have poor SFS 
model performance and, in fact, have instabilities in simulations with peri- 
odic boundary conditions [5] . This result highlights the importance of cor- 
rectly interpreting norm regularizations. To illustrate this, we compute 
256^ DNS and 64^ models runs for their case 3a: decaying Taylor-Green 
vortex with f = 2.5 • 10"'^ (time step, dt = 3 ■ 10~^ for all simulations). In 
Figs. [3] and SI we display results for the models employing a filter width 
a — 27r/16 (results with a = 27r/24 are closer to that of the DNS). Fig. 
[3^ shows the time evolution of the total kinetic energy in the super-filter 
scales. 



where ka = 27r/a. Incorrectly applying the norm, E — |v|^/2, to 
Clark— a and LANS— a leads to the interpretation that they deviate signif- 
icantly from the DNS due to slower energy decay. Correctly applying the 
norm, Ea ~ v ■ v/2, shows that the models are, in fact, doing quite 
well as can be seen in the plot of the dissipation rate of energy, e (Fig. 
[SJj). The qualitative properties of Lagrangian models is that they have 
the same invariants as the primitive equations, albeit in a different norm. 
Since it is these invariants which may very well infiuence the dynamics, as 
for example shown in |12) where Kolmogorov spectra are present at inter- 
mediate times at large scale in ideal (Euler) three-dimensional fiuid fiows, 
one must compare the behavior of one set of invariants to the other set 
in the Lagrangian modeling formulation. Similarly, the norm spectra 
for Clark— a and LANS— a in Fig. appear under-dissipative at high 
wave numbers, while the correct norm spectra are closer to the DNS. 
We also conduct Clark— a and LANS— a computations for case 3b of [S]: 
V = 1/3000, a = 27r/32 and 384^ grid points (Fig. gb). We find no signs 
of instability. The particular numerical expression of the models we used 
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Figure 2: (a) Compensated energy spectra (27r/a, vertical dashed line) for 
Navier-Stokes (solid black), LANS— a (blue dash-dotted), Clark— a (green 
dashed), and Leray— a (red dotted), (b) Normalized structure function 
scaling exponent ^p/fs versus order p. Clark— a is the best approxima- 
tion for the superfilter-scale spectrum, whereas high-order intermittency 
properties are best reproduced by LANS— a 



8 



0.0 



2 4 6 8 10 

t 

0.10 1 .- J 




0.00 r , , , I , , , I , , , I , , , I , , , 1 

2 4 6 8 10 

t 



Figure 3: CaseSa, Re w 400. Line styles are as in Fig. [2^ for correct norms 
with the addition of norm energy for Clark— a (pink dash-triple-dotted) 
and for LANS— a (cyan long-dashed), (a) Filtered kinetic energy, tke. (b) 
Dissipation, e = dE/dt, versus time, t. 



is given in Eqs. ^ and ©. These differ from Eqs. (8) and (9) in [5] 
and it is known that discretized pseudospectral operators depend on their 
algorithmic form (e.g., Vi/'^ 7^ 2ijj'V'4' [131 HI)- This is hkely the source 
of their observed instabihty; however a recently discovered deficiency in 
high-order low-storage Runge-Kutta schemes [B] may also impact numer- 
ical implementation of the models. For this reason, our calculations are 
made with second-order Runge-Kutta in time. It is worth pointing out 
that Clark— a is nearly identical to the Rational Large Eddy Simulation 
(RLES) model [TS] which has been shown to be unstable for non-periodic 
boundaries [SJ [32] ■ Both our work and that of jS] have employed periodic 
boundary conditions only, but extra precautions must be taken for stability 
in the non-periodic case (Sj [32] . 



3 MHD: circulation and outlook for LES 

In MHD, the circulation properties are quite different since small-scale 
circulation is broken by the Lorentz force j x b, with j = V x b the current, 
b being the induction. This force acts as source (sink) of circulation, F, 
as opposed to the insufficient modeling of F in the Leray— a and Clark— a 
models. This can be seen in Kelvin's theorem, 

4 /v-dr= /j xb-dr. (11) 



dt dt jQ jQ 

This may prevent the formation of rigid bodies even while conserving all 
the correct physical properties of the flow. 
The LES equations for MHD are given by 

9tv -fa)xv==jxb — VH -I- vV^w — V • f 

(9tb = V X (v X b) + ryV^b - V • f\ (12) 

where 77 is the magnetic diffusivity, H = P + |bp/2 the modified pressure, 
f is the Reynolds SFS stress tensor, 

= - ViVj - (bibj - hibj) , (13) 

and is the electromotive-force SFS stress tensor, 

= hvj - hvj - {v,bj - Vibj) . (14) 

Note that the extension of eddy viscosity to eddy resistivity employs the 
usual form for f involving only the filtered velocity while the expression 
for similarly only involves the filtered magnetic field [52]. Meanwhile, 
Eqs. ([rs]) and ([T4| make it explicitly clear that interactions between the 
two fields at subGlter scales must be taken into account. 
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Figure 4: (a) CaseSa, Re ~ 400: energy spectra after peak of dissipation, 
t G [4.5, 5.5]. Line styles are as in Fig. [H Vertical dashed line is 2TT/a. (b) 
Case3b, Re « 3000: dissipation, e, versus time, t for Clark— a (light green) 
and for LANS-a (dark blue). 
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Another problem with extending eddy-viscosity concepts to MHD is 
that they can be related to a known power law of the energy spectrum . 
This is inappropriate for MHD as neither kinetic nor magnetic energy is 
a conserved quantity and the general expression of the energy spectrum 
is not known at this time [SI |M1 [H [HI ISn]- Additionally, MHD has 
been shown to have nonlocal interactions between large and small (su- 
perfilter and subfilter) scales [1] (e.g., in the Batchelor viscous-inductive 
regime, u ^ rj, where energy is transferred directly from viscous-scale, l^, 
eddies to small-scale magnetic field, I [5]). This complex interaction 

is a challenge in general for MHD-LES, but may be an advantage for the 
Lagrangian-averaged approach as energy exchanges with sub-a scales may 
disrupt rigid body formation. Some limited case MHD LES include the 
cross-helicity model [48] assuming alignment between the fields and the 
low magnetic Re LES [Ml 35] ■ Extensions of spectral models to MHD 
based on two-point closure formulations of the dynamical equations pro- 
posed recently look promising in the analysis of turbulent flows and of 
the dynamo mechanism [3]. Approximate deconvolution models for MHD 
|36| also appear promising, but have yet to be tested on non-laminar flows. 
However, there are many regimes of MHD dependent on the ratios between 
the various conserved quantities and i' /rj: there has yet to be demonstrated 
a generally applicable LES for MHD. 

3.1 LAMHD— a and absence of rigid bodies 

The Lagrangian-averaged MHD a-model (LAMHD-a) gSlllTJIlT] is given 
by, where the velocity if filtered as before and b = (1 — a;^V^)~^b: 

9tV -|-a;xv=jxb — Vtt + i^V^v 
Sfb ^ V X (v X b) + ryV^b 
V-v = V-v = V-b = V-b = 0. (15) 

LAMHD— a may be written as a MHD-LES, Eqs. (HI]), for the case aM = Oi 
[23l[25], which we study here. The model preserves the ideal quadratic in- 
variants of MHD (in the li\ norm) as well as Alfven's theorem for frozen- 
in field lines [28]. Physically, it supports Alfven waves at all scales while 
slowing and hyper-diffusively damping waves with wavelengths. A, smaller 
than a [5S]. In examinations of its SFS physical properties LAMHD— a 
(dashed lines) displays neither positive-exponent power-law scaling nor 
superfilter-scale spectral contamination (see Fig. [5]). Under similar condi- 
tions LANS— a (not shown) displays these signs of rigid body formation. 
Further examinations with larger filters and higher Re were unable to un- 
ravel any sign that rigid bodies form for LAMHD— a [25^ . 
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Figure 5: Kinetic (a) and magnetic (b) energy spectra {2TT/a vertical 
dashed line). The largest scales are affected by differences from the MHD 
DNS (solid lines) in initial conditions. LAMHD— a (dashed lines) exhibits 
neither the positive power-law nor the superfilter-scale spectral contamina- 
tion associated with high Re LANS— a. 
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3.2 LAMHD-a as a SFS model 



Since LAMHD— a did not display any signs of rigid body spectral contam- 
ination, we test it as a SFS model for large kinetic and magnetic Reynolds 
numbers, « 3300. A DNS of MHD is computed at a resolution of 1024^. 
The initial conditions for v and b are a superposition of ABC modes ]W\ 
with random phases and wavenumbers k G [1)4]. No external forcing is 
applied and the total energy is allowed to freely decay. LAMHD— a is com- 
puted for identical conditions at a resolution of = 168^ with the same 
Re and a filter size a = 6 A = 27r/28. As a base- level comparison we also 
compute an under-resolved (or no-model) solution of the MHD equations at 
A^^ = 168^. Time evolution of the total energies and enstrophy are shown 
in Fig. |6l In comparison with under-resolving MHD, LAMHD— a shows 
errors of approximately the same magnitude in these global quantities. 
Comparisons of energy spectra (Fig. [7]), however, show an improvement 
in predictive quality for LAMHD— a, especially for the magnetic energy 
spectrum. As turbulence develops, energy begins to pile up at small-scales 
and deplete at intermediate scales for 168^ MHD. LAMHD— a improves 
the prediction of superfilter-scale spectra compared to no SFS model. 

Cross-sections of shown in Fig. |5]at t — 8.4 indicate that LAMHD- 
a finds sharper and better defined, more intermittent current sheets than 
the under-resolved run compared to the DNS. 

4 Summary 

Incompressible LANS— a, while it performed well at moderate Reynolds 
number is limited as a high Re SFS model. Due to its strong suppres- 
sion of spectrally local interactions at subfilter-scales, and consistent with 
its conservation of small-scale circulation, LANS— a develops rigid bod- 
ies which contaminate the superfilter-scale energy spectrum. In contrast, 
Clark— a and Leray— a, neither of which conserve small-scale circulation do 
not develop energy-spectrum contamination from rigid bodies. LANS— a, 
however, best matches the intermittcncy properties of Navier-Stokes fluid 
turbulence. 

In MHD, a mechanism for local small-scale transfer is the interaction 
of small-scale Alfven waves. As LAMHD— a supports Alfven waves at all 
scales while slowing and hyperdiffusively damping those with wavelength 
A < a, it more gently suppresses SFS local interactions than LANS— a. 
This together with the greater nonlocality in MHD and the Lorentz-force 
source of small-scale circulation, inhibits the formation of rigid bodies in 
LAMHD— q;. It appears to retain the good intermittency properties of 
LANS— a without its poor spectral properties. For this reason, we find 
LAMHD— a to be a viable model at high Re in 3D. As LAMHD— a has been 
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Figure 6: Time evolution of total (solid), magnetic (dashed), and ki- 
netic (dotted) energies (a) and total dissipation, e, (b) for 1024'^ DNS 
(black/solid), 168^ LAMHD-a (blue/dash-dotted), and 168^ no-model 
(red/dotted). LAMHD— a provides no improvement in prediction of global 
quantities over an under-resolved DNS. 
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Figure 7: Compensated kinetic (a) and magnetic (b) energy spectra aver- 
aged over t G [8,8.8]. Line styles are as in Fig. |6l (27r/a vertical dashed 
line). Energy piles up at small scales in the no- model approach (under- 
resolved DNS) and LAMHD— a is seen to provide improved predictions of 
the superfilter-scale spectra, especially for the magnetic field. 
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Figure 8: Cross sections of square current, j^, at t = 8.4 for no-model 168"^ 
(left), 168^ LAMHD-a (center), and 1024^ DNS (right) at full resolution. 
LAMHD— a provides a much better capturing of the intermittent current 
sheets than the under-resolved solution. 

previously found to reproduce the difficult to model properties of MHD at 
high Re in 2D [45^ and moderate Re in 3D [44] , we believe it will prove to 
be a generally applicable MHD LES, in many instances in geophysics and 
astrophysics where magnetic fields are known to be important dynamically. 
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